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PHASED MODELS FOR EVALUATING THE 


PERFORMABILITY OF COMPUTING SYSTEMS 

by 

L. T. Wu and J. F. Meyer 


The University of Michigan 
Ann Arbor, MI 48109 


Abstract - On-line control applications of fault-tolerant com- 
puters often require the computers to execute different sets of 
computational tasks during different phases of a control process. 

To evaluate the system's "ability to perform", a phase-by-phase 
modeling technique is introduced. Intraphase processes are allowed 
to differ from phase to phase, where the probabilities of interphase 
state transitions (which occur at the time of a phase change) are 
specified by interphase transition matrices. Based on constraints 
imposed on the intraphase and interphase transition probabilities, 
various iterative solution methods are developed for calculating 
system performability . 


I . INTRODUCTION 

During recent years, the use of probabilistic models as a basis 
for evaluating the performance and reliability of computing systems 
has become increasingly widespread. Typically, the models employed 
are Markov processes (e.g., II]) or queueing models (e.g., [2]) which 
can often be analyzed in terms of imbedded Markov processes. However, 
it is usually assumed that the underlying process is "time-homogeneous" 
in the sense that state transition probabilities are invariant with 
time. (By "state" here we mean the state of the "total system", i.e., 
the state of the computing system and its environment.) Although this 
assumption of time-homogeneity is appropriate for certain applications, 
there are many situations where the user's demands on the computing 
system can change appreciably during different phases of its utiliza- 
tion. This is particularly true for real-time control applications in 
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which the computing system is required to execute different sets of 
computational tasks during different phases of a control process. 

One approach to dealing with a time-varying environment is to 
decompose the system's utilization period into consecutive time peri- 
ods {usually referred to as a decomposition of the system's "mission” 
into "phases"; see [3]-[5J, for example). Demands on the system are 
then allowed to vary from phase to phase; within a given phase, 
however, they are assumed to be time-invariant. This permits intra- 
phase behaviors to be evaluated in terms of conventional time- 
homogeneous models, but raises the interesting question of how the 
intraphase results are combined. This is the essential question 
addressed in investigations of "phased mission" reliability evalua- 
tion methods (e.g., [3] -[5]) where the problem has been constrained 

as follows. It is assumed, first, that a "success criterion" {formu- 
lated, say, by a "structure function"; see [5] for example) can be 
established for each phase, where the cirterion is independent of 
what occurs during other phases. It is required further that success- 
ful performance of the system be identified with success during all 
phases, that is, the system performs successfully if and only if, 
for each phase, the corresponding success criterion is satisfied 
throughout that phase. 

Although the above constraints are reasonable for certain types 
of systems, they exclude systems where successful performance involves 
nontrivial interaction among the phases of the mission. In more 
exact terms, it has been shown (see [6] , Theorem 6) that such 
"structure-based" formulations of success are possible if and only if 
the phases are "functionally independent" in a precisely defined 
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manner. What we v;ish to do, therefore, is to examine the utility 
of "phased models" in a less restricted context. 

In addition to removing the above constraints , we extend the 
domain of application to include evaluations of computing system 
"performability" [7], [8]. (Although performability concepts will 
be introduced as needed in the presentation that follows, some prior 
familiarity with this background may improve the reader's perspective 
of what is being accomplished.) Finally, unlike the models used in 
phased mission reliability evaluation, we permit the state sets of 
the intraphase models to differ from phase to phase. Thus, the 
modeling of a particular phase can be tailored not only to the 
computational demands of that phase but also to the relevant pro- 
perties of the total system which influence performance during that 
phase. 


II. PHASED MODELS 

Generally, in modeling the performability of a computing system 
C in some specified computational environment E (see [7], [8]), the 
most detailed view of the total system S = (C,E) is represented by 
a stochastic process X g referred to as the base model of 5 . X g is 
defined over a time interval T called the utilization period and each 
random variable X^ (teT) takes values in a state space Q, i.e., with 
respect to a common "description space" ft, X t :ft->-Q. In general, a 
state qeQ represents a particular status of both the computer C and 
its environment E. Moreover, the computer coordinate of q may include 
both the structural state of C and the internal state of the struc- 


ture. 
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An instance of the base model's behavior corresponding to an | 

outcome w is a state trajectory (or "sample path") where f 

u : T-j-Q with u (t) = x. (u>) , for all teT. The collection jl 

td U) u 

tj 

u = 

i' 

is referred to as the trajectory space of S . At a higher, less fj 

f! 

detailed level of description, the user's view of total system behavior 

tj 

is modeled by a random variable Y^, called the performance of S. Y„ ij 

\D " _ O 

takes values in a set A of accomplishment levels where 't is assumed 
that Xg is refined enough to support Yg, i.e., there exists a function 
Yg :U-)-A, called the capability function of S , such that, for all wefi, - 

Y S = Y S ^ ' 

Finally, the performability of S is taken to be the probability dis- 
tribution function of the performance variable Yg or, in case A is 
discrete, the probability mass function p g where, for all aeA 
Pg(a) = the probability that S performs at level a. 

To generalize the notion of a "phased mission" in the context of 
performability modeling, let us suppose that the utilization period 
T is the continuous interval T = [0,h]. Suppose further that T is 
divided into a finite number of consecutive phases (time intervals) 

f 2 ~ t 2 ~ [ti 7 12 ^, * • • ' T m ~ f tm- 1 ' ^"m^ where 0 — t q ^ ^ ••• 

< t m = h. During phase T^, we assume that the system can be modeled 
in the manner described earlier for the entire period T, i.e., by a 
(continuous time) stochastic process 

X k = {X k |teT k > 

where each random variable X^ takes values in the phase k state space 
0 k (X fc *. fl->-Q k ) . X is referred to as the intraphase process (of phase 
k) and, combining these processes, we obtain the process 


m k m . 

X = U X = U {X?|teT.}. 
s k=l k=i c Jc 

On examining x g we see that it is similar to a base model except 

that, for each time instant t^, l£k£m-l,the state of the system is 

k k+1 

represented by two random variables X. and X, whose values, 

k . . z k 

respectively, are the final state of the k phase and the initial 
state of the k+1 phase (see Figure 1) . However, if we consider 
an augmented utilization period 

T = T u ( t£ | k=i ,2, , m-1 ) 

(where t^ can be interpreted as the initial time of pahse k+1) , 
then Xg can be expressed as 
X g = {X t |t£T> 

where 

if t = 0 

if te(t k-i'V u) 

i f t = 

If, further, we regard the state space of X g as the union 
m 

Q = U Q k 
k=l K 

then Xg is a base model in the sense defined above. When X^ is so 
constructed from intraphase processes, we will refer to it as a phased 
base model. 

Let us suppose now that the base model X g of a performability 
model is phased and that X g "supports" the capability function y g in 
the sense that the end-of -pahse "samples" of a state trajectory u 
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uniquely determine the accomplishment level y s (u). More precisely! 
a phased model X g supports y g if 


u(t.) = u’(t.), for all k, implies y q (u) = y q (u'). 


( 2 ) 


Given that X g supports y g , the performability model can then be sim- 
plified as follows. The simplified base model is taken to be the 
imbedded discrete-time process 

X s = {z k |k=0,lf . . . ,m} 


where Zq = Xq and, for each k>l, 
z k = <_ • 


(3) 


Since Z Q is required only for the initial state distribution, the 
trajectory space U of 5Tg can be effectively regarded as the product 


space 


U = Q 1 X Q 2 x ’ * 


* Q. 


m 


(where is the state space of phase k) . The corresponding simpli- 
fication of y is the capability function 
y s :U A 

where if u(t^) = q^, for k=l, 2 , . . . ,m, then 

V9a/g 2 V = V u) • 

Then, by (2), it follows that, for all aeA, 


p s (a) = Pr(Yg 1 (a)) = Pr^ta)) 


(4) 


and hence the performability model (X ,y s ) can be used to evaluate 
the performability of S. We will thus refer to (Xg,y s ) as being 
equivalent to the model {X g ,y s ) . 

Although the concept of "support" (2) might appear to be somewhat 
restrictive, this is not the case when we look at what is typically 




done in reliability modeling. Given a traditional single-phase 
reliability model, th ® system reliability can often be determined 
by ‘sampling the state of the system at the end of its utilization 
period. Such single phase equivalents (or multiple phase equiva- 
lents in the case of phased models) exist whenever traditional reli- 
ability modeling assumptions are made with regard to the intraphase 
processes. 

To illustrate this point, consider a continuous time Markov 
model of a TMR (Triple-Modular-Redundancy ) system with a perfect 
voter where the simplex system has failure rate X, i.e., the base 
model Xg = {X t |tsT} is represented by the graph 

1 9 

2 A 3A 

3 A 2A . 

If the utilization period is T = [t^,t^] and the accomplishment set 
is A = {ag,a^} (where a^ = success and a^ = failure) then a^ is 
accomplished if and only if the system is in state 1 or 2 throughout 
T. Thus the capability function is 

1 if u(t)e{l,2}, for all teT 

Y s (u) - 

0 otherwise 

and accordingly the performability at a Q (i.e., the reliability) is 

P s (a Q ) = PrfYg 1 ^)) = Pr({u|y s (u) = a Q }) . 

Since state 3 is absorbing, it follows that Pr({ujy s (u) ^ a^ and 
u (t^) e(l, 2} )) = 0, and hence 

P s (a Q ) ~ Pr ({u| Y s (u) = a Q }) + Pr({u|y s (u) ^ a Q and uft^ e{l, 2} }) 
= Pr ( { u j u ( t^ ) f: { 1 , 2 } } ) 

= Pr(X t e{l,2}) . 


Accordingly, the equivalent base model is a pair of random variables 
X s = {Zp-Zj.) 

describing the state of the original model at the beginning and the 
end of the utilization period, i.e., Z n = X, and Z, = X. . Tl.e 
corresponding equivalent capability function is the structure function 
7 S :Q A where 

a Q if qe{l, 2} , 

7 s (q) = 

' a^ otherwise. 

III. PROBABILITY COMPUTATION OF CARTESIAN TRAJECTORY SETS 

If (X c ,y c ) is a perforniobiiity model of S, then the performability 
of S is determined by the probabilities of the trajectory sets Yg (a)?D 
(see ( 4 )) where, for each aeA, Yg^ ( a ) the set of all state trajec- 
tories of Xg corresponding to accomplishment level a. Generally, the 
evaluation of Prfyl^a)) requires a detailed knowledge of how intra- 
phase processes cooperate to accomplish level a, i.e. , a thorough 
understanding of their "functional dependencies" (see [6]). The dif- 
ficulties are further aggravated by statistical dependencies between 
phases. However, we have found that when a trajectory set VSJ is 
Cartesian in the sense that, for every phase k, there exists R^pQ^ 
such that v = R-, x r„ x . . . x r , then Pr(V) can be determined iter- 
atively using matrix multiplications. Moreover, given this ability to 
compute the probabilities of Cartesian sets, the probabilities of more 
general sets can be determined by decomposing them into disjoint 
unions of Cartesian components. (The latter is taken care of auto- 
matically by algorithms which determine Yg'*’ (a) ; see [8])* Hence, the 
problem reduces to that of computing the probabilities of Cartesian 
trajectory sets. 


If X g is the phased model from which X g is derived, for each 
phase k, let = X*. t be the initial state of the k intraphase pro- 
cess and let n^ be the number of states in Q^. Then, for a Cartesian 
trajectory set V = R 1 x R 2 x ... x r , the conditional intraphase 
transition matrix of the k phase is the n, x n,_ matrix P,, where, 

‘ — rn ™ K K V / K 

for all i,j€Q k > 

P V,k {i,j) = Pr(2 k =j l Y k =i ' Z k-l eH k-l' •••' z ie R i ) • 

In other words, P . (i,j) is the initial -to-final state transition 

V / K 

til * * • 

probability of the k, intraphase process, conditioned by the first 

k-1 components of V. Similarly, for all but the first phase, the 

conditional interphase transition matrix is the n^_^ x n^ matrix 

H v , k where, for all ieQ k-1 and jeQ k , 

H V,k^'^ = Pr ^ Y k“^l Z }.-l =i ' Z k-2 eR k-2' * ’ * ' z i eR i^ • 

i. 

In other words, H . (i,j) is the probability that the k phase 

V / K 

th 

initiates in state j given that the final state of the k-1 phase 
is i, conditioned by the first k-2 components of V. For consistency, 
we let Hy ^ be the n-^ x n^ identi.ty matrix. Finally for each phase, 
the characteristic matrix of the k th phase is the n^ x matrix 

G V,k where 


! 1 if i=j and ieR^ 

0 otherwise. 

In terms of the above matrices, we are able to establish the following 
matrix formula for computing the probability of a Cartesian trajectory 
set V. Given X g , let 1(0) denote its initial state distribution, i.e., 
1(0) = [p 1 ,p 2 , . . . ,p n J where p i = Pr(Z Q =i) = Pr(xJ=i); and let F k 
denote the x 1 column matrix with "1" in each entry. Then by 
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induction on k f it can be established that 
Theorem 1 : 

If V = R, x ... x r x Q. , x ... x Q then 
l k k+1 m 

Pr(V) = I(0)-[J i H V/J ,.P v>i .G V;1 ].F k . 

Proof : 

For k»l, 

I(0) * H V,1’ P V,1 = I(0) ’ P V,1 = [a lf . . . ,a jr . . . ,a n ^3 

where 

n. 


3 j = Pr(Z 0 =i) -Pr ( 2i=j | 2 0 =i) 


n. 


= > i>r (Z Q =i,Z 1 =j) = Pr (Z^=j ) . 

i— 1 

Multiplied by ^ and F^, 

I ( 0 >- H v , r p v , r G v , i- F i 

= J PrfZ^j) = PrUj^) 
j SR 1 


= Pr(Z 1 eR 1 ,Z 1 eQ 2 /.. .,Z eQ m ) 

= Pr (V) . 

Suppose that the formula holds for k<m, then 
k+1 

I (OH n H„ n -p„ „-G rT J -F 


1=1 


v,i v i a k+i 


H V' S '' P V> st ‘' Gv ’^ ,H V,k+l ’ P V,k+X‘ G V,k+l * F k+l 


A 1 * H V, k+1 - P V, k+1 - G V, k+1 * F k+l 
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where 

^7 = [b , t ... ) b ] 

3 ”x 

and 

Pr ( Z k = j ' Z k-l eP ‘k~l' * • * ' z 1 eR 1 > if 

0 otherwise, 
by applying the equation for k. 

When we iteratively compute the matrix product, beginning from 
the left, the first two terms become 

A 2 ~ A l’ H V,k+l ~ l c i • • • • /c., . . . , c ] 

J n k+l 

where 



n. 


'j ~ ^ hi’^k+l* 1 '^ 

iL Pr (z k = 1 ' z k-l £R k-l' • * * ,z i eR i ) * 


Pr ^ Y k+l ^ Z k ^^k-l^k-l' ' * * ,Z l eR i) 


« l 


ieR ^ ( Y k+l~ j ' V 1 ' Z k-l £R k-l " • • ' Z 1 £R 1 ) 
k 

Pr (Y k+l =:j ,Z k £R k' * * ' ' z i ei V * 

The next partial product is the result of multiplying A 2 by the 
transition matrix which yields; 

A„ = A-* P,_ 


3 A 2‘ F V,k+l = . ,d . , . . . ,d ] 

J k+1 


where 


n. 


d. = 


k+l 

l 

i=l 


C i* P V,k+l (l, 3) 
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n 


‘k+l 


i=l PrtY k+l a * i,z k eR k r * * * ' z i eI V 

Pr ^ Z k+l~^ l Y k +l” i ' Z k e R k ' • • • /Z 1 eR i ) 


n 


l k+l 

ill Pr (Z k+l =j ' Y k+l =i ' Z k eR k' . . . ,Z 1 er jl ) 

= Pr (Z k+l =j ' z k eR k ' ' * ' • 

The product is completed by multiplying a 3 by the characteristic 
matrl * ° f the k+lth Ph3Se *" d «» —i»9 vector P k+1 , that is , 

H v,Ji ,p v,r 0 v, 2 J - p k+ i 
A 3 ’ G v,k+i' p k+i 

i£R Pr (Z k+l~ 1 /Z k eR k' * * ■ ,z i eR i) 

K+X 

Pr ^ Z k+l eR k+l' Z k eR k' ’ • * • z i eR i^ 

1 Rl# ’ * * ' Z k+l ER k+l' Z k+2 eQ k+2' * ‘ * ,Z m eC V 

— Pi"(R,x...xr xq x vn l 

i k+l u k+2 * ' ’ Q m J • 

Accordingly, the equation holds for all k < m , whioh completes ^ 
proof of Theorem 1. 

Corollary : 

For any Cartesian set v = H xR x... XR , 

A ^ m 

, m 

Pr(V) = 1(0). [ n H .p r i 

1=1 V '£ v, z ' G v, a ^ ’ F m* 


C5) 


illustrate thus method, consider a system with three identical 

of the subsystems is dedicated to different 01 

different computational tasks. 
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However, during the second phast ~ [t^t^], the system is recon- 
figured into a TMR configuration. The system is capable of degraded 
performance which occurs when (i) at least one subsystem has failed 
during phase 1 and at least two subsystems are functional through. u; 
phase 2, or f ii ) no failures occur during phase 1 and the system 
functions as a simplex system at the end of phase 2. Suppose that 
each of M^, M 2 and M^ fail permanently with a constant failure rate 
X and failure characteristics of the subsystems are statistically 
independent. Then the probabilistic nature? of phase 1 and phase 2 
can be represented, respectively, by finite-state time-homogeneous 
Markov processes with transition graphs as illustrated in Figure 2. 

Based on the above description of the system, three accomplish- 
ment levels can be established, i.e., A 55 {a Q ,a^,a } where = no 
failure, a^ = degraded performance and a 2 = failure. When expressed 
in terms of the state trajectories 

Yg- l (a 0 ) = {(1,1) , (1,2)} = {l}x{1.2} 

Y^ 1 (a 1 ) = {2,3,4}x{l,2}u{l}*{3} 

7^ 1 (a 2 ) « {2 ,3,4,5}x{3,4}[j{l}x{4}u{5}x{l,2} 


Then, solving the intraphase probabilities which, in this case, are 
the same for all Cartesian sets V, 


V,1 


r 

0 

0 

0 

0 


2 

r s 


0 

0 

0 


2 

r s 


0 

2 


0 

0 


2 

r s 


0 

2 


3rs 2 +s 3 


2rs+s‘ 

2rs+s‘ 

2rs+s' 1 

1 


r = e 
s = 1 - r 


and 
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V,2 


P 

0 

0 

0 


3p 2 g 


p 

o 

o 


3pq' 

2pq 

P 

0 


q 

q 

l 


p = e 
q = 1 - p 


Suppose that the initial state distribution of X g is 1(0) = [1,0, 0,0,0]. 
If at time t 1 , the system is in state 5 with respect to the phase 1 
model (i.e., at least two subsystems have failed) then, depending on 
the exact number of subsystems failed, the state of the system with 
respect to phase 2 model is either 3 or 4 . 

By applying the definition of interphase transition probability, 

Pr(two failures before t^) 

H V , 2 ^ 5 ' 3 ^ Pr(two or three failures before ~t^) 

-Ut -t ) 

- 3e P _ 

-A(t r tJ C 1 
1 + 2e 1 U 


and 


H V,2 {5 ' 4) = 


1 - e 


“X (t-tg) 


-Mt^tg) C 2 


1 + 2e 

Transitions from states other than 5 happen to be deterministic, and 
thus we obtain the following interphase transition matrix 


V 


1 

0 

0 

0 

0 


0 

1 

1 

1 

0 


0 

0 

0 

0 

c- 


0 

0 

0 

0 

c. 


Ky by definition, is the 5 *5 identity matrix. 
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Using these intraphase and interphase matrices, the probability 
of V = y^( 3 q) can be computed using equation 5, i.e., since 



1 0 0 0 0 - 
0 0 0 0 0 


~1 0 0 0 

G v,l ~ 

0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 

and G Vt2 = 

0 10 0 
0 0 0 0 
0 0 0 0 


we have 

Pr(V) = 1(0) * H v,l’ P V,l' G V,l' H V,2* :E> V,2' G V,2‘ P 2 
= r 3 p 3 + 3r 3 p 2 q . 

Similarly, 

Pr(y^ 1 (a 1 )) = 3r 2 sp 2 + 3r 3 pq 2 

and 

Pr ( a 2 ) ) = s 3 +3rs 2 +3pqr 2 s+3qr 2 s+r 3 q 3 . 

Although equation 5 provides us with a general formula for computing 

the probability of a Cartesian set, its disadvantages derive from 

the fact that H , and P . matrices may be difficult to obtain in 

V / JC V / JC 

practical applications. In particular, these matrices will generally 
depend on V as well as X g and, moreover, will generally depend on the 
history of X g prior to phase k. However, the latter objections dis- 
appear when the transition probabilities are "memoryless". More pre- 
cisely, let the (unconditional) intraphase transition matrix of the 

Vi 

k phase to be the n k x n k matrix F k where, for all i,jeQ k , 

P k (i,j) = Pr(Z k =j|Y k =i) , 

1.U 

i.e., the probability that the k intraphase process ends up in state 
j given that it initiates in state i. Similarly, let the (unconditional) 
interphase transition matrix be the n k _^ x n k matrix H k where, for all 
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iEQ k-l and ^ eQ k' 

H k (i,j) = p r (^ k =j I z j c -. 1 = 3 i) ' 

■f* In 

i.e., the probability that the k intraphase process initiates in 
• ■ th 

state 3 gxven that the k -1 intraphase process ends up in state i. 
Then the intraphase transitions of X„ are memoryless for V at phase k 

b 

if 


P = P 
V,k k 

Similarly, the interphase transitions of X g are memoryless for V at 
phase, k if 

H v,k = V 

Accordingly, when transitions are memoryless through phase k, by the 
definitions and Theorem 1 we obtain 
Theorem 2 : 

If V = r ^ x P’ 2 >< * • • xR k*Qk+l X ’ ' " *®m and ^- n ^ ra Ph ase an< ^ inter- 

phase transitions of X g are memoryless for V through phase k, then 

Pr(V> = KOWJ . 


Corollary : 

For any Cartesian set V, if the intraphase and interphase tran- 
sitions of X are memory less for V for all phases, then 
s 

m 

Pr(V) = ICOM^ . 

When V is a Cartesian set and , for H - 1,2,..., k-1, then the 

intraphase and interphase transitions of Xg are memoryless for V 
through phase k. Moreover, G v ^ is an identity matrix for l - 1,2,..., 
k-1. Accordingly, applying Theorem 2, we obtain the following formula 
for' the probability of the trajectory set V = . . . x Qj < _x ><R k x Q ] c +7 * * * * x Q m 



-17- 


which, alternatively/ is the probability of the event "Z^eR k ". 
Theorem 3 : 

If V = Q^ x . . . x Q k _^ x R k x Q k _^ 1 * * * * then 

k 

Pr(V) = I(0).[ n h £ .p a j .G v ^.F k . 

Z ” 1 


When Theorem 3 is specialized to singleton sets R k = {i}, where 

ieQ k , it permits us to compute the probability of the event "Z^i 1 

More generally, if we denote the probability distribution 

of the random variable Z k by the n k ~dimensional vector I (k) = 

[p ,P 2 ,..., P ] where p i = Pr(Z k =i), then 

k 

Theorem 4 : 


k 

l (k) = I(0M n H P ‘P P 3. 

z=i * * 


Proof : 

By Theorem 3, 

Pr(Z k =i) = Pr (Q 1 x . . . x Q k _ 1 x{i}xQ k+1 x, . . xq^) 


k 



where, by the nature of V, G v k *F k =E^, i.e., the n k * 1 column matrix 
with "1" on the i n entry and "0" elsewhere. Thus, Pr(Z k =i) is equal 
to the i n entry of the n k ~dimensional row vector 
k 

I(0H n H . -P.] . 

£,=1 * * 

By Theorems 2-4, when certain intraphase and interphase transi- 
tions are memoryless for V, the probability of a Cartesian set V is 
relatively easily obtained. However, such results may still be 
difficult to use due to the fact that, even though the transitions 
are msmoryless for V, they may not be memoryless for other Cartesian 
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sets, Accordingly, we have sought to identify stronger conditions 
under which the formulas will hold for all Cartesian trajectory sets. 
First, by extending previous definitions, the intraphase (interphase) 
transitions of X g are memoryless at phase k if they are memoryless 
for all Cartesian sets V at phase k; the intraphase (interphase) transi 
tions of X„ are memory less if they are memoryless at all phases. The 
advantage of memoryless transitions are obvious, for by their defini- 
tion and the corollary to Theorem 2, we have 
Theorem 5 : 

If X is a phased model and the intraphase and interphase transi- 
tions of Xg are memoryless then, for all cartesian sets V, 

m 


Pr (V) = I(0M n 
£=1 


W G V,£ ] ' F m‘ 


( 6 ) 


Moreover, we find that the memoryless property is relatively easy to 
characterize, that is, we are able to show the following characteris- 
tic conditions for the memoryless property. It is important to note 
that the conditions do not involve any specific Cartesian sets. 

T heorem 6 : 

(1) The intraphase transitions of X g are memoryless at phase 
k if and only if, for all i^eQ^ (£=1,2, . . . ,k-l) , 

fr(z k =j|Y k =i,z k _ 1 =i k-1 ,...,Z 1 «i 1 ) =P k (i,j). (7) 

(2) The interphase transitions of X g are memory less at phase k 
if and only if, for all i^eQ^ (£“1, 2, . . . ,k~2) , 

Pr(Y k =j|Z k-l =i,Z k-2 =i k-2 Z l =i l ) ” H k (i ' 3) ' (8) 

Proof ; 

Suppose P„ , is memoryless for all Cartesian sets V = R. *R *...XR 

V / X Z ] 

By taking R^ to be the singleton set {i^}, £ = l,2,...,k-l, 

P V,k (i '^ = Pr ^k = ^ I Y k =1,Z k-l =:L k-l' * ' * /Z l =i l ) 


= P k (i# j ) • 
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Now, suppose that, for all i^eQ^ U = 1, 2, . . . ,k-l) , 

Pr ( Z k = j I Y k = 1 ' Z k-l = 1 k-l ' * * * /Z l = 1 l ) 

= Pr(2 k =j|Y k =i). 

Then, for any Cartesian set V - R,xR„x...xR , 

l 2 m 

P V,k* i,: ^ = Pr ^ Z k = ^ Y k = 1 ,Z k-l eR k-l / * - ’ ' z l eI V 


i eR I i eR CPr(Z k j l Y k“ i/Z k~l =i k-l' 
1 l eR l' * * * /:L k-l eR k-l 


* Pr(Y k =i ' Z k-l =i k-l 


Pr (Y^-i , f ■ • •/Z^gR^T 

Thus, by the assumption, P v k (i,j) is equal to 


l Pr (Z k -j |Y k -i) -Pr {Y k -i,z k _ 1 -i k ^ 1 , . 

:L 1 £R 1 / • v* ,J 'k-l eR k-l 

p r( y k =i,Z k _i R k _ 1 '**wZ 1 R ] _) 

Factoring out the term Pr (Z k =j | Y k =i) , we have 
P V,k (i/ ^ = Pr ( z k =j l^ k =i) • 1 = P k (i,j) 


• ' z l“ 1 l ) 

• • , Z i~i ^ ) I 


' Z l =i l 


) 


which completes the proof for part (1) of the theorem. Part (2) is 
proven in a like manner. 

Finally, when the transitions of X g are memoryless, the equiva- 
lent model Xg is a time- varying Markov chain. This eyn be demonstrated 
as follows. 

By Theorem 6 and the definition of memoryiess transition, we have 


Pr(Z k 1 k ^ Y k ^ ' Z k-1 1 
= Pr{Z k =i kl Y k^ } 

= Pr(Z k =i k |Y k =j,Z k _ 1 


k-1 ' ' * * 


=i k-l> • 



Accordingly, for each k<m and i^Q^ U - l,2,...,k). 
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Pr(Z k ^^k-l^k-l" * * /Z l‘' 1 l ) 

= £ Pr(Z k =1 J V^^k-l^k-l' * * * ' 


IjeQ 


Pr( Y j l z k-rvi z i =i i ) 


= ,L pr( ^ i ki^'Vf i k-i )>pr( v^ z k-fVi ) 


jeQ 


ieo Pr ^ 2 k _1 k' Y k~^ ^ z k-l i k-l ) 

“ Pr (Z k =1 kl Z k-l =1 k-l ) ‘ 

Hence, X g = {z k | k= 0 , 1 ,. . . ,m} satisfies the Markov properties. 

Moreover, the transition probabilities of X g associated with 
phase k can be expressed as a matrix. 

P k = f ? k (i,j)] 


where 

P k (i,j) = Pr(Z k -j 12^=1) 

- X pr( v j lv h ' z k- 1 = i) - pr! v h l z k-i =i) 

hsQ k 

= l Pr(Z k =j|V k -h) -Pr{Y =h|z. ,=i) 
heQ k x k k xi 

= l P, (h, i) ■ H (i,h) . 
h£Q k K k 

Accordingly, in terms of matrix multiplication, 

- v p * 

and equation 6 can be represented in a more convenient form: 

m _ 

Pr(V) = KOU ^ P r G v>t ).P m . 


- 21 - 


IV. SUMMARY 

It has been shown that the concept of a "phased mission' uan be 
extended to performability models via the notion of a "phased" base 
model X g . Under reasonable conditions, X g yields an equivalent per- 
formability model (Xg,Yg) and, as demonstrated by the results of the 
paper, the intraphase and interphase probabilities of X g suffice, to 
determine the probabilistic nature of X g . In particular, it has been 
shown that, for any trajectory set V of X g with a "Cartesian" struc- 
ture, the probability of V can be computed as a product of matrices 
(Theorem 1). In general, each matrix depends on X g and V but, as 
established in subsequent results (Theorems 2-5) , the formulations 
may be simplified when certain phases of x g are "memoryless for V." 
Finally, it has been demonstrated (Theorem 6) that transitions which 
are memoryless for all Cartesian sets V are characterized by a 
"Markovian property" relative to preceding end-of-phase observations 
of the phased base model X g . 
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